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ABSTRACT 

The distributions of Secchi depths (i.e. water 
transparency) with simultaneously measured standard oceano- 
graphic parameters on file at the National Oceanographic 
Data Center are surveyed on a global basis. An inventory of 
many of the oceanographic parameters is given for all Marsden 
squares. Correlation coefficients between Secchi depths 
and eleven other parameters are also tabulated. Linear 
regression equations for some twenty-one selected ocean areas 
relating Secchi depth and the other parameters are presented, 
and in some cases plotted. No simple and consistent relations 
between Secchi depth and other parameters are evident; how- 
ever, several trends are noted. Porel color and oxygen 
measurements show trends toward an inverse proportionality 
with Secchi depths while bottom depth data indicate a 
possible direct proportionality. 
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I. INTRODUCTION 



A . GENERAL 

The optical properties of sea water govern a great portion 
of the nature of the oceans of the world. Light penetration 
below the sea surface is of interest in a large number of 
areas, since it is due to this characteristic that animal 
and plant are allowed to exist in the oceans. Strangely 
enough, studies of the optical nature of sea water seem to 
have received less emphasis than the extensive investigations 
by chemists and biologists. In years to come, optical ocean- 
ography will play a larger role in our biological, physical 
and military studies of the oceans. 

Arsen'yev and Voytov (1968), in their studies of Bering 
Sea water, noted the possibility of an inverse relationship 
betvreen the level of plankton biomass and sea water trans- 
parency. Uda (1963) found that transparency (i.e. Secchi 
depth) and dissolved oxygen measurements could be helpful in 
investigating productivity levels. Kampa and Bowden (1957) 
have shown that vertical migrations of deep scattering layers 
depend directly on light penetration below the sea surface. 
Voitov and Dement' yeva (1971) analyzed the relationship 
between transparency and the plankton concentration in 
Antarctic and Indian waters. They concluded that the depth 
to which the Secchi disc can be seen underwater is a function 
of plankton content and suspended terrigenous material, the 
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distribution of which is linked with hydrodynamics. Ryther 
(1969) has estimated that over half of the worlds fish supply 
is produced in no more than 0.1# of the ocean. This fact, 
coupled with the necessity of light energy for growth of 
plants in the sea will make future research in optical 
oceanography of greater significance to the biologist. 

An important part in the understanding of the oceans on 
a global basis lies in the realm of optical oceanography. 

For example, Jerlov (1968) has proposed a scheme of classi- 
fication of ocean waters of the world into several oceanic 
and coastal optical groups by irradiance characteristics. 

This was an important step in optically grouping similar 
waters over the world. 

Militarily, sea water optics has always been important. 

It is vital for a submarine commander to consider the water 
transparency, since his vessel' s hull is continually susceptible 
to visual detection by patrol aircraft and surface ships. 

Boden, Kampa and Snodgrass (i 960 ) and Kampa (1970) have re- 

* 

ported an apparent correlation between extinction coefficients 
and thermocline depths. These depths at which significant 
thermal changes take place are of considerable importance to 
buoyancy control in a submersible. A slight change in the 
light conditions below the surface can also affect visual 
mine detection and diver operations as well as determine the 
success or failure of underwater photography. 

Serving as the primary source of energy for organisms in 
the oceans, light supports photosynthesis and permits life 



9 



in the sea. Through the study of light penetration * man 
may someday achieve greater control of the immense potentials 
of the oceanic harvests. Murphy (1959) describes this effect 
of water transparency on the albacore fishing industry. 

For purposes of discussion in this paper, "transparency" 



refers to the average of the depths at which the Secchi disc 
disappears upon lowering and reappears upon raising. "Color," 
or "Forel color," is the color code obtained by the comparison 
of the color of the Secchi disc when lowered to a depth of 
one meter with a standard color scale such as the Forel-Ule 
scale. Beam attenuation coefficient (a) is a measure of 
the attenuation of a collimated light beam through a fixed 
path length. Irradiance is a measure of the light flux per unit 
area striking a horizontal subsurface plane. Vertical ex- 
tinction coefficient (K.), or diffuse attenuation coefficient, 



is a function of wavelength and is defined as follows: 




where 



K 



X 



vertical extinction coefficient [meters 



I 



irradiance at the sea surface 



o 



I 



z 



irradiance at depth = Z meters 



B. PURPOSE 



The purpose of this investigation is to analyze Secchi 
disc measurements and, on a selected area basis, determine, 



if possible, correlations between these and other simul- 
taneously measured standard oceanographic parameters. This 
was accomplished using the Naval Postgraduate School IBM 
System /360 Model 67 computer and the Biomedical System 
library available. Global data was obtained from the 
National Oceanographic Data Center, Washington D. C. 

C. THE SECCHI DISC 

Over the years, the Secchi disc has become a widely used 
tool by which oceanographers and marine biologists have 
measured water transparency. Commander Cialdi and Professor 
Secchi (1968) made first mention of its use after observa- 
tions at sea in 1865 . Secchi analyzed, how the disc visibility 
changed as a function of its color, the solar altitude, sea 
surface reflections, ship’s shadow, sky clearness, water 
color, and the height of eye above the water surface. En- 
countering near perfect sea conditions, he observed an increase 
in disc visibility with increased solar altitude and disc 
whiteness. 

Poole and Atkins (1929) measured submarine illumination 
in the English Channel and developed a w T idely used empirical 
formula relating Secchi depth (Z ) and vertical extinction 

o 

coefficient (K) : 



K = 1 . 7/Z [rrf 1 ], 

O 
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This is equivalent to 



* 

where Z is measured in meters, 
s 

saying that the illumination at the Secchi depth is about 
16# that penetrating the water's surface. Gall (1949) 
independently confirmed this in separate readings ranging 
from 16# to 2 5 55 » Bumpus and Clarke (1947) in their studies 
on oceanic, coastal and harbor water also provided confirma- 
tion of this relationship. Murphy (1959) j for example, 
estimated the ecological effect of water transparency on 
albacore catches using this equation. Holmes (1970) has 
mentioned that although the Secchi depth measurement is not 
precise and is dependent upon many variables, it is still 
useful as a practical measure. For instance, biologists 
interested in primary production have come to use three 
times Secchi depth as an estimate of the depth of the 
euphotic zone. Holmes (1970) proposed the Poole-Atklns 
constant of 1.7 should be reduced to 1.44 for turbid water 
calculations. Russell and Clarke (1944a, 1944b), in a 
series of cruises, measured transparency in a number of 
oceanic and coastal areas. Among items they noted to affect 
transparency were bottom depth, zooplankton and phytoplankton 
content, water mass movement and harbor runoff. Atkins, 



* Henceforth K is written without the subscript X. It 
is to be understood that the dominant wavelength in the case 
of K obtained from Secchi data is determined by a number of 
factors, including the spectral attenuation of sea water, 
the spectral distribution of incident solar radiation, and 
the spectral acuity of the observer's eye. 
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Jenkins and Warren (195^) examined suspended matter in sea 
water and its seasonal effect on Secchi depth readings. 

They oberserved an inverse relation between the amount of 
phytoplankton and the visual range of the disc. Graham 
(1966) compiled and analyzed Secchi depths, vertical extinc- 
tion coefficients and Porel color data from several cruises. 
Multiple regression formulas were obtained expressing the 
vertical extinction coefficient as a function of color and 
the reciprocal of Secchi depth. Graham noted the excellent 
agreement between vertical extinction coefficient and Secchi 
depth but added caution must be heeded when extrapolating 
from one area to another. Visser ( 1 9 6 7 ) examined the effect 
of wave action on disc visibility for different disc sizes 
and concluded a fair relationship between observed mean 
Secchi depth in meters and the corresponding x> e :rcentage 
yellow in the Porel color according to the formula: 



100 



.26Y + 1.9, 



where Y represents the percent yellow of the \v'ater. 

Tyler ( 1 9 6 8 ) discussed the theory of how the Secchi depth 
can be used to calculate the sum of the vertical extinction 
coefficient (K) and the beam attenuation coefficient (a) 
utilizing the formula: 

3. 78 K _ Z s 

Tzrr-KT - z 1Q% > 
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where Z^o% re P resents the depth in meters at which 10% of 
surface illumination was measured. 

Holmes (1970) followed Tyler’s recommendations and found 
no change in disc visibility as a result of varying disc size. 
He also concluded that using a viewer below the sea surface 
made no difference in relatively turbid water. Holmes sug- 
gested that the pulsating subsurface light field created as 
a result of wave action, when combined with vertical disc 
movement, was responsible for compensating for the loss of 
contrast caused by glare and refraction at the surface. 

Holmes, too, was interested in the correlation between 
Secchi depth and the thickness of the euphotic zone. Visser 
( 1969 ) as well confirmed the biologists' euphotic zone rule- 
of-thumb and, in addition, independently confirmed the 
Poole-Atkins relationship. 
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II. DATA PROCESSING 



A. TREATMENT OF DATA 

Oceanographic data from some 86,258 stations were received 
on magnetic tape from the National Oceanographic Data Center. 
The information included all data recorded on the station 
A-sheets at sea including all chemistry. Limited computer 
storage available necessitated screening data to eliminate 
storing measurements not to be used in this study. The 
original data, however, have been preserved and remain at 
the Naval Postgraduate School for future work. Station data 
selected for analysis were : 

Sec chi Depth 
Latitude 
Longitude 
Marsden Square 
Water Depth 
Forel Color 
Cloud Cover 
Month of year 
Sample Chemistry: 

Temperature 

Salinity 

Sigma- t 

Oxygen 

Phosphate 

Phosphorus 

Nitrite 

Nitrate 

Silicate 

Sample pH 



15 



Data obtained but not utilized in this study are: 

Cruise Identity Code 
Station Identity Code 
Country Identity Code 
Date/Hour of Sample 

Number of Observations at this Station 

Wave Height/Direction 

Wind Direction/Speed 

Barometric Pressure 

Weather Information 

Air Temperature 

Sound Velocity 

Due to computer storage limitations, only sea surface 
data were analyzed. Following disc storage of the data, it 
was necessary to scan by Marsden Square and tabulate all 
measurements available for analysis. The results of this 
survey are included as Table I and serve as a reference as 
to what data are available worldwide. Data were extremely 
sparse beyond Marsden Square 571; therefore, they were not 
analyzed. It is seen from Figure 1 and Table I that the 
great majority of Secchi disc measurements are concentrated 
along the shipping routes of the world. Worthy of mention 
is the lack of open ocean data and the abundance of 
measurements in Japanese coastal waters. 

Upon completion of the initial data survey, it was 
necessary to select a number of areas of interest for analysis. 
The selection process was based on the following three consid- 
erations. First, attention was given to areas of high data 
density, since a more representative analysis results from 
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a larger data sample ; second, consideration was given to 
areas of greater oceanographic interest such as upwelling 
regions and major current systems, and an attempt was made 
to insure there was a complete cross section of all water 
types; and finally, emphasis was given areas of some military 
interest. Using these criteria, 44 areas were selected, 
and correlation coefficients and cross tabulation plots were 
obtained for each using the computer program BIOMED 02D.. 

To aid in the analysis of data distribution, all data 
w T ere further broken down by one-degree Marsden sub-squares. 
This tabulation was used to analyze as much data as possible 
in as small an area as was reasonable. In this way unneces- 
sary averaging of data from varying water types was avoided. 

At times, especially in areas near Japan, it was necessary 
to select segments of high data Marsden squares for analysis, 
since attempting to analyze data for all squares would have 
exceeded computer storage limitations. For example, nearly 
700,000 bytes were required to analyze area 7 (Shikoku Basin). 
Likewise, when data were sparse, several Marsden squares were 
combined to form an analysis area. Table III is a summary 
of correlations obtained for the 44 areas chosen, and Figure 
2 shows parameters from several areas plotted for comparison. 

After linear plots and correlation coefficients were 
obtained for all areas, 21 regions were selected for further 
analysis using the BIOMED 02R least squares approximation 
linear regression program. Figures 3 through 44 show some 
of the more successful results of these analyses and are 
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discussed in Section III of this report. The selection of 
the 21 areas for regression analysis was based on visual 
inspection of the linear plots for some evidence of a consis- 
tent relationship between the various parameters and Secchi 
depth. Daniel and Wood (1971) list several curve types which 
were used for comparison as were several hand plotted examples 
of logarithmic and linear relationships. Figure 3 shows a 
BIOMED 02D graphical output w ith the corresponding regression 
curve from BIOMED 02R, plotted. Appendix B is a sample output 
of the package program BIOMED 02R and Appendix A is an example 
of the Fortran program utilized to screen data for the BIOMED 
programs . 

Chemistry data at times were too sparse to allow analysis. 
Table I shows that although temperature and salinity data 

are usually abundant, such is not the case for measurements 

# 

of oxygen, phosphate, nitrate, nitrite, and silicate. In 
light of this, wherever data were plentiful enough, a special 
effort was made to attempt correlations. The graphs from 
all areas were compared, and special notice was made of any 
consistencies. All graphs were analyzed as to special area 
considerations. For example, Figure 14 shows a plot of 
salinity and Secchi depth. The wide variation of salinity 
measurements is presumably due to fresh ’water runoff into 
the area. Figure 26 shows water depth plotted against Secchi 
depth in which transparency is seen to increase with depth. 
These are examples of special area factors which were 
considered in the analysis of the graphs and correlation 
coefficients . 
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B. BIOMEDICAL COMPUTER SYSTEM PROGRAMS 



1. General 

The Biomedical (BIOMED) computer programs were 
developed at the University of California at Los Angeles 
(Dixon, 1971) and designed for simplicity and convenience. 
Little previous computer experience is necessary as the pro- 
grams are in "package" form and require just a few control 
and system cards to activate. 

2 . BIOMED 02D (Correlation with Transgeneration) 

This program is intended to be basic in data sorting 
and analysis. Pairs of data, or cases, are selected for 
inclusion in analysis if they satisfy a specified Boolean 
statement. This permits exclusion of undesirable or blank 
data and allows creation of data subsets. Mathematical data 
manipulations, or transgenerations, are also a feature of 
the program. Output from this program includes: sums, 

means, cross product deviations, standard deviations, 
variance-covariance matrices, Correlation matrices, and one 
page cross-tabulation plots of any two variables. The data 
density code utilized in BIOMED 02D graphical outputs is 
translated in Table VI. 

3 . BIOMED 02R (Stepwise Regression) 

This program computes a sequence of multiple linear 
regression equations in a stepwise manner. At each step, 
the variable which makes the greatest reduction in the error 
sum-of-squares is added to the regression equation. The 
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variable added is the variable which has the highest partial 
correlation with the dependent variable partialed on the 
variables that have already been added, and equivalently, it 
is the variable which, if it were added, would have the 
highest F value defined below. Variables can also be forced 
into the equation. 

For each step in the stepwise regression procedure, 
the following are computed: 



df 


= 


n - q 


a dd 


= 


E (X k(J - X)(X k<J - X d ) 
k=l 


b dd 


= 


i=n ? 

SS = Z (y, - Y, T 

i=l 1 1 


MS 


= 


SS/df 


rdf 


= 


q 


RSS 


= 


a dd “ b dd 


RMS 


= 


RSS/ rdf 


F 


rr 


RMS/MS 



where : 



n 




total number of samples 


q 


= 


number of variables involved 


df 


= 


residual degrees of freedom 


SS 




sum of squares, residual 


MS 




residual mean square 


rdf 


= 


regression degrees of freedom 


RSS 




regression sum of squares 
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RMS 



regression mean square 

the value of the independent variable 

the mean value of the independent variable X^ 
th 

the i value of the dependent variable 
t h 

the i value of the least squares approximation 
of the dependent variable 



X, 



kd 



X 



d 



Y. 

1 



As indicated by the above equations, a higher P value 



shows that a higher degree of confidence may be placed in 
the presented regression equation. Generally speaking the 
F value means little until compared with other similarly 
determined values. The user may also indicate specific 
numerical values of P which act as limits to allow or prevent 
the variables from entering the regression equation. Appen- 
dix A shows the first two steps of a BIOMED 02R regression 
analysis output. Boolean case selection available in BIOMED 
02D is not a feature in BIOMED 02R. Therefore, it was neces- 
sary externally to screen rough data prior to calling upon 
BIOMED 02R for regression analysis (see Appendix A). At 
each step in the procedure the multiple correlation coeffi- 
cient R serves as an indication of how well the regression 
equation fits the data. With R = ± 1.0, a perfect correlation 
is indicated. A widely used measure of data fit is the square 
of the multiple correlation coefficient (R ). Another option 
which can be included in the BIOMED 02R output are plots of 
the residuals (computed Secchi depth minus actual Secchi 
depth) against the independent variables in the equation. 
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Analysis of these plots aids in the determination of closer 
data fits. 

A wide selection of transgenerations is available 
to alter data sets (or subsets) to more desirable forms 
before inclusion in the analysis. 
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III. DISCUSSION OP RESULTS 



A. DATA SURVEY 

Table I shows the data breakdown for Marsden Squares 1 
through 571. Frederick (1970) surveyed Secchi disc data and 
reported that a great deal of the available data were from 
Japanese waters. This is apparent in Table I which shows 
large concentrations of data in Marsden Squares 9^-97 and 
129-132. 

Figure 1 and Table I, however, indicate mid-oceanic data 
to be sparse. Since a study such as the present one must 
consider all variables affecting transparency measurements, 
mid-oceanic data would be valuable since the additional 
dynamic coastal processes tend to complicate matters. If 
a consistent relationship between an oceanic parameter and 
Secchi depth is to be found, it would probably be most appar- 
ent away from seasonal coastal effects such as fresh water 
runoff and upwelling. Thus, the lack of mid-oceanic data 
reported to date is unfortunate. 

Secchi measurements are also rare in the polar Marsden 
squares, and for that reason data are tabulated only through 
Marsden Square 571. Table IV presents parameter means by 
area. 

B. CORRELATION COEFFICIENTS AND LINEAR REGRESSION ANALYSIS 

Tabulated results of BIOMED 02D correlation analyses for 

the 4A areas appear in Table III and for some areas are plotted 
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in Figure 2. Table V is a summary of the regression analysis 
results (BIOMED 02R). Area delineations were made according 
to Figure 1 and Table II. Figures 3 through 44 are examples 
of the BIOMED 02D graphical output. Alpha-numeric data 
density codes are translated in Table VI. A correlation 
coefficient of ± 1.0 indicates perfect correlations, while 
a value of 0 signifies no correlation at all (Table III). 
Correlations between Secchi depths and each parameter will 
be discussed separately. 

1 . Color 

It is to be noticed, first of all, that Forel color 
code seems consistently to correlate better than other 
parameters, as indicated by higher coefficients in Figure 2 
and Table III. It is seen that these coefficients are 
usually negative, signifying lower transparencies for higher 
color code. A similar pattern was reported by Visser (1967). 
Figures 3 through 9 show this inverse trend for seven areas, 
where Forel color is plotted against Secchi depth. An 
inverse, possibly exponential, trend is unmistakable in 
Figures 3 through 8; however in Figure 9 for area 35 (Eastern 
coast of the United States) it is evident this trend is 
interrupted near colors 7 through 11. Figure 1 shows area 
35 to be an area of high fresh water runoff as well as an 
area of high industrial activity, which accounts for possibly 
high terrigenous and chemical influences on water color 
measurements. Several attempts at regression equations were 
made to relate color code and Secchi depth. These equations 
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are presented by ocean area in Table V. Figure 3 shows two 
regression equations for area 7 graphically. Notice that the 
first approximation (Curve A) achieved a multiple regression 
coefficient of 0.65, while for the second there is a coeffi- 
cient of 0.72. Since data density in area 7 (Shikoku Basin) 
was high, numerous other regressions were attempted as shown 
in Table V. One of the most successful of these resulted 
from noticing the apparent exponential nature of the color- 
Secchi depth graphs. The results were the equations: 

Step 1 

In (Z ) = 3.73 - 0.21} (Color) (R = 0.70) 

Step 2 

In (Z ) = 5.05 - 0.17 (Color) 
s 

- 0.3 (Oxygen[ml/l] ) (R = 0 . 76 ) 

with the multiple regression coefficients listed. Table III 
shows that higher multiple regression coefficients did result 
from this study; however, none were as consistently high as 
those for color code. 

2 . Sea Surface Temperature 

Examination of Table III yields no obvious trends in 
correlations between sea surface temperature and Secchi depth 
however, several interesting cases were encountered. Figures 
10 through 12 give several of these. Figure 10 shows tempera 
ture and Secchi depth plotted for area 35 (Eastern coast of 
United States). Notice the large temperature range (25°C) 
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as compared to Figure 11 for area 12 (Western Pacific Ocean). 
Figure 1 shows area 35 to be a coastal area, affected by the 
coastal influences discussed previously for Figure 9* Area 
12, on the other hand, is seen to be more oceanic with a 
resulting smaller temperature range (10°C). A regression 
line is plotted in Figure 10 on the temperature data from 
area 35, resulting in a somewhat low multiple regression 
coefficient of 0.36. 

3. Salinity 

Nearly all the areas analyzed were coastal areas, 
subject to localized salinity variations caused by fresh 
water runoff. This effect can be seen in Figures 13 through 
19. In each case, the lower salinity water is associated with 
decreased transparency. This result is not surprising 
considering the high levels of terrigenous suspensoids in 
river waters. Figures 20 and 21 are included for comparison, 
since they represent somewhat oceanic areas (9 and 12), away 
from land drainage. The observations are reasonable in that 
the high salinity variation is not present, and the data 
are therefore more closely grouped. 

Figure 22 shows the salinity and Secchi disc data 
for area 30 (Peruvian coast). Note that both disc visibility 
and salinity data variances are small. 

Plotted in Figure 23 are Mediterranean Sea data 
which display the very high salinities characteristic of 
that area. No significant data correlations resulted for 
the area. 
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Superimposed on Figure 18 is the regression line 
attempted for area 36 (Labrador Sea). This line depicts 
the equation: 

Z g = 2.18 + .01 (Salinity [°/ 00 ]) 2 

which results in a multiple regression coefficient equal to 
O. 85 . Due to the extremely local nature of the fresh water 
runoff, this result would, of course, be of little 
significance in other coastal areas. 

4 . Bottom Depth 

Bottom depth is closely related to water transparency. 
As deeper water is encountered away from the coast, three 
depth-related influences play contributing roles. First, 
shallow water transparency is determined to a great extent 
by bottom agitation. The absence of stirred-up bottom sedi- 
ments in deep oceanic areas cause water transparency to in- 
crease. Secondly, the coastal effect of river runoff causes 
decreased water transparency since these rivers carry large 
quantities of both suspended and dissolved materials. These 
suspensoids and solutes may then be carried hundreds of miles 
by coastal currents and cause great increases in coastal water 
turbidity. Finally, coastal regions which lie at the eastern 
extremities of oceans tend to be upwelling areas. In these 
areas, nutrients from deeper water are carried to the surface 
and a reduction of Secchi depth results with enhanced 
phytoplankton growth. 
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These influences of water depth on Secchi depth can 



be seen in Figure 2 and Table III, where nearly always a 
positive correlation coefficient results between bottom 
depth and water transparency. This indicates an increase in 
bottom depth is accompanied by an increase in Secchi depth, 
an anticipated outcome in view of the factors discussed 
above. The results of this line of reasoning can be seen 
graphically in Figures 24 through 28, where data from five 
areas show this direct relationship. Figure 27 shows a 
regression equation plotted for area 13 (Bay of Bengal), 
for which a high multiple regression coefficient was attained. 
Although this trend between bottom depth and Secchi depth 
was prevalent, it was not always the case as can be seen in 
Figure 28, where an upwelling area (Equadorian coast) is 
depicted with bottom depth and water transparency apparently 
unrelated . 

5 . Density 

No real trend is evident in the sigma-t results 
presented in Figure 2 and Table III. Figure 29 through 31 
are graphs of three coastal areas where the data distribution 
looks similar to salinity plots given in Figures 13 through 
19* Fresh water runoff is postulated as the reason for this 
in all cases. 

6 . Oxygen 

Figure 2 shows a negative correlation nearly always 
resulted when measured values of oxygen level and Secchi 
depths were compared; however, no consistent relationship 
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was developed. Figures 32 and 33 are graphs of oxygen level 
and Secchi depth. A regression equation for area 34 (Sar- 
gasso Sea) is plotted on Figure 29 yielding a somewhat high 
multiple correlation coefficient of 0.74. Figure 30 shows 
the high variability of oxygen data from coastal area 7 
(Shikoku Basin). Time variations in oxygen level can be 
expected due to river runoff and photosynthesis-respiration 
cycles . 

7 . Total Phosphorus 

Since phosphorus is basic to the growth of phyto- 
plankton in the sea, seasonal and geographic variations of 
phosphorus level in the coastal areas analyzed can be expected. 
As phosphorus is returned to the sea upon the death of plants 
or organisms, it would possibly serve to increase light 
scattering and thereby lessen penetration depths. This 
resultant decreasing effect on Secchi depth measurements is 
evident in Table III and Figure 2, where nearly all correla- 
tion coefficients are negative. No consistent correlation 
level was observed since, in most areas, data were too sparse 
to permit analysis. Figure 34 shows total phosphorus data 
plotted. 

8 . Silicate 

Silicate level can be expected to be generally higher 
in coastal waters due to the effect of land drainage. Although 
silicate measurements are somewhat sparse, negative correla- 
tions did prevail in the areas analyzed as shown in Table III 
and Figure 2 . 
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Figure 35 shows silicate data for area 28 (Northwestern 
coast of United States), in the vicinity of the Columbia River 
discharge. There is an apparent exponential decrease in 
transparency level as the silicate level increases. Super- 
imposed on Figure 35 is a plot of the equation which was 
obtained from regression analysis. 

9. Nitrate, Nitrite and Phosphate 

Figures 36 through 38 are graphs of nitrate, nitrite, 
and phosphate data, respectively. Without exception, data 
were too sparse to permit any form of accurate correlation 
or regression analysis. This lack of data is shown in the 
breakdown of data, Table I. 

10 . Cloud Cover 

No correlation was noted between the amount of cloud 
cover and Secchi depth for any of the 44 areas analyzed. 

Figure 39 shows cloud cover data for area 26 (Gulf of Alaska). 

11. Latitude and Longitude 

Analysis of Secchi depth data as it varied with 
latitude and longitude yielded no significant correlations 
or regression results; however, at times the graphical re- 
sults were of interest. An example is Figure 40 (Black Sea), 
where latitude and Secchi depth data are plotted. The in- 
fluence of river discharge in the northern Black Sea is clearly 
apparent as water transparency decreases significantly. On 
the other hand. Figures 4l and 42 show that for the two up- 
welling areas (both South American coastal areas), transparency 
is somewhat independent of latitude changes. Figure 43 
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shows the variation of water transparency with longitude 
for area 13 on the western extremity of the Indian Ocean. 

As longitude is increased, greater influence from coastal 
effects, such as bottom stirring, upwelling, and land drain, 
are to be anticipated. The data distribution in Figure 4 3 
was therefore, not surprising. In contrast, Figure 44 is 
a typical data distribution for an area where water trans- 
parency is independent of longitude; accordingly. Figure 1 
shows area 24 (Aleutian Islands) to be away from these 
coastal influences. 

C . MONTHLY AVERAGES 

Monthly Secchi depth averages were plotted in Figure 45 
for areas 13 (Bay of Bengal)and 35 (U.S. East Coast). Figure 
46 shows Secchi depth averages for upwelling areas 28 
(Northwestern coast of United States) and 29 (Peruvian coast). 
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IV. CONCLUSIONS 



Secchi depth readings are influenced by many sea water 
parameters; among them are sea surface temperature, salinity, 
density, water depth, color, oxygen, phosphorus level and 
silicate level. The exact nature of these relationships 
cannot be determined on the basis of the present study; 
hov/ever, analysis of correlation coefficients yields general 
trends. Of particular consistency are the high correlation 
coefficients relating Secchi depth and Forel-Ule color code. 
Several of the parameters analyzed yielded consistently 
negative or positive correlation coefficients. 

The data survey reveals that water chemistry records in 
most areas of the world are too sparse for linear regression 
analysis techniques. Among these are nitrate, nitrite, and 
phosphate measurements. 

Coastal areas under the influence of land drainage 
present wide ranges of temperature, salinity, density and 
chemistry readings. 
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V. SUGGESTIONS FOR FURTHER RESEARCH 



Work should continue to provide better worldwide coverage 
of Secchi disc measurements and of the other parameters 
normally sampled. 

The scope of this project was limited to sea surface 
measurement analysis because of computer storage requirements. 
Further analysis should continue to analyze subsurface 
parameters and their possible correlation to Secchi disc 
measurements . 

The Secchi disc and procedures for its use should all 
be standardized to eliminate possible variations in future 
data due to differing disc dimensions and measurement 
practices . 

As mid-ocean data become available, other studies such 
as this should be made, since correlations determined free 
from coastal influences may yield results which can be 
simply and accurately extrapolated to other open ocean areas. 

There exists in use a package program similar to, but 
superior to, the BIOMED programs. Known as the "Statistical 
Package for the Social Sciences" or "SPSS" (Nie, Bent, and 
Hull, 1970), this fully integrated system allows for subfile 
creation, missing values and conditional transformations. 

An analysis using SPSS would be superior and much simpler to 
any BIOMED-based study. 



33 




3'l 



DATA DISTRIBUTION BY MARSDEN SQUARE 
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TABLE I. (cont’d) MARSDEN SQUARE DATA DISTRIBUTION TABLE 
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TABLE I. (cont’d) MARSDEN SQUARE DATA DISTRIBUTION TABLE 
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TABLE II 



AREA DELINEATION TABLE 

AREA ENCOMPASSED 

Marsden Square 201 

Marsden Square 200 

Marsden Square 199 

Marsden Square 167 

Marsden Square 165 

Latitude: 3^ ° — 38° 

Longitude: 123°-127 

Latitude: 30°-35° 

Longitude: 135°-l40 

Also included area 
bounded by: 

Latitude: 35°-36° 

Longitude: 139°-l40 

Marsden Square 130 

Marsden Square 129 

Marsden Square 97 

Marsden Square 96 

Marsden Square 95 

Marsden Squares 63 

Marsden Square 62 

Marsden Square 6l 

Marsden Square 60 

Marsden Square 59 



(N) 

° (E) 
(N) 

° (E) 

(N) 

° (E) 
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AREA NO. 



AREA ENCOMPASSED 
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Marsden Square 24 
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20 


Marsden Square 22 
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Marsden Squares 360-361 
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Longitude: l60°-l80° (W) 
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Latitude: 50°-55° (S) 




Longitude: 60°-70° (W) 
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Marsden Square 43 
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AREA ENCOMPASSED 



AREA NO. 
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Marsden Square 115 
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Marsden Square 186 


37 


Marsden Square 216 


38 


Marsden Squares 214,215 




Latitude : 


60°-65° (N) 
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Marsden Squares 327-328 
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TABLE III. PARAMETER CORRELATION COEFFICIENT SUMMARY 
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TABLE III. PARAMETER CORRELATION COEFFICIENT SUMMARY (Cont’d) 
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TABLE III. PARAMETER CORRELATION COEFFICIENT SUMMARY .( Cont ’ d ) 
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Figure 2A Correlation Coefficient Graph 
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Figure 2B Correlation Coefficient Graph 
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Figure 2C Correlation Coefficient Graph 
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Figure 2D Correlation Coefficient Graph 
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TABLE IV (Cont'd) 

Parameter mean summary 
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TABLE V LINEAR REGRESSION ANALYSIS RESULTS 
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TABLE V (Cant'd) LINEAR ESGPZSSION ANALYSIS RESULTS 
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TABLE V (Cont’d) LINEAR REGRESSION ANALYSIS RESULTS 
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TABLE VI 

GRAPH DATA DENSITY CODE 



This following table shows the symbol representations used 
in plotting frequencies of data in BIOMED 02D graphical output. 
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Figure 15. 
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Figure 17 
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Figure 18. 
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Figure 19. 
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Figure 20 
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Figure 2h 
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Figure 26. 
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Figure 29 
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Figure 30. 
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Figure 32 
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Figure 33 
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Figure 35* 
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Figure 33 
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Figure liO. 
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Figure lj2 
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Figure h3« 
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Figure hii 
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SAMPLE BIOMED 02R OUTPUT( cont 1 d) 
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